Analytic solutions in non-linear massive gravity 
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We study spherically symmetric solutions in a covariant massive gravity model, which is a candi- 
date for a ghost-free non-linear completion of the Fierz-Pauli theory. There is a branch of solutions 
that exhibits the Vainshtein mechanism, recovering General Relativity below a Vainshtein radius 
given by {rgni?Y^^, where m is the graviton mass and rg is the Schwarzschild radius of a matter 
source. Another branch of exact solutions exists, corresponding to Schwarzschild-de Sitter space- 
times where the curvature scale of de Sitter space is proportional to the mass squared of the graviton. 
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Introduction: It is a fundamental question whether 
there exists a consistent covariant theory for massive 
gravity, where the graviton acquires a mass and leads to 
a large distance modification of General Relativity (GR). 
The quest for massive gravity dates back to the work by 
Fierz and Pauli (FP) in 1939 [1]. They considered a mass 
term for linear gravitational perturbations, which explic- 
itly breaks the gauge invariance of GR. As a result, there 
exist five degrees of freedom in the graviton, instead of 
the two found in GR. There have been intensive studies 
on what happens going beyond the linearised theory. In 
1972, Boulware and Deser (BD) found that, at the non- 
linear level, there appears a sixth mode in the graviton 
that becomes a ghost in the FP model [2] . This problem 
was reexamined using the effective theory approach [3], 
where additional (Stiickelberg) fields were introduced to 
restore the gauge invariance, and whose scalar part rep- 
resents the helicity-0 mode of the graviton. In the FP 
model, the scalar acquires non-linear interaction terms 
that contain more than two time derivatives, signaling 
the existence of the ghost. 

The Stiickelberg approach also sheds light on the other 
puzzle in the FP gravity: if one linearises the system, the 
solutions in the FP theory do not reduce to GR solutions 
in the massless limit. This is known as the van Dam, Velt- 
man, Zakharov (vDVZ) discontinuity [4, 5]. However, 
in this massless limit the scalar mode becomes strongly 
coupled and one cannot linearise the system. Therefore, 
due to strong coupling, the scalar interaction is shielded 
and GR can be recovered. This is known as the Vain- 
shtein mechanism [6]. The strong coupling scale in the 
FP model is identified as A5 = (m^Mp/)^/^ where Mpi is 
the Planck scale and m is the graviton mass. This scale is 
tightly connected with the non-linear interactions of the 
scalar mode that contain more than two time derivatives. 
In the decoupling limit, where m — >■ and Mpi — > 00, 
while the strong coupling scale A5 is kept fixed, one ob- 
tains an effective theory for the scalar mode, where it is 
possible to study the consistency of the theory in more 
detail. 

Until recently, it was believed that there is no consis- 
tent way to extend the FP model [7, 8] to get a ghost 



free model at all orders. A breakthrough came with a 
5D braneworld model known as Dvali-Gabadadze-Porrati 
(DGP) model [9]. In this model there appears a contin- 
uous tower of massive gravitons from a four dimensional 
perspective. The non-linear interactions of the helicity- 
mode of massive gravitons contain no more than two 
derivatives, which is crucial to avoid the BD ghost. Due 
to this fact, the strong coupling scale in this theory is 
given by A3 — {m?MpiY/^ instead of A5, where m = 
and Tc is a cross-over scale between 5D and 4D gravity 
[10, 11]. Further studies have considered more general 
non-linear interactions which contain no more than two 
derivatives. In 4D, only a finite number of terms satisfy 
this condition; these are dubbed Galileon terms because 
of a symmetry under field transformations of the form 
9^j7r — > dfj_-K + [12]. Ref. [13] constructed the exten- 
sion of the FP theory that gives the Galileon terms in 
the decoupling limit, by choosing the correct parameters 
in the Lagrangian up to quintic order in perturbations. 
Ref. [14] proposed a covariant non-linear action that au- 
tomatically ensures this property to all orders, which we 
will discuss below. 

A remaining crucial question is whether this property, 
holding in the decoupling limit, is sufficient to ensure the 
absence of the BD ghost or not. In Ref. [14], it was shown 
that there is no BD ghost in the decoupling limit to all 
orders in perturbation theory, but only up to and includ- 
ing quartic order away from this limit. However, it is 
very hard to show the absence of the BD ghost at all or- 
ders if one starts from Minkowski and studies non-linear 
interactions perturbatively. Therefore, it is important 
to obtain non-perturbative background solutions in this 
theory, and study fluctuations around them. Moreover, 
it is interesting to find solutions in this covariant non- 
linear theory, that can describe features of the observed 
universe. These are the topics of the present work. 

Covariant non-linear massive gravity: We first 
construct the action for generalised FP model based on 
Ref. [13, 14]. We define the tensor H^i, as a covarianti- 
zation of the metric perturbations: 



V^lidf^rducb^. (1) 
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The Stiickelberg fields = (a;" — tt") transform as 
scalars, while ijap corresponds to a non-dynamical back- 
ground metric that is needed to define the potential, 
which is assumed to be the Minkowski metric. The co- 
variant tensor iJ^^ can then be expanded as 

= V-2m-' (2) 

and under the coordinate transformation — -|- 
TT^ transforms as 

7r'^^7r^+^^. (3) 

Before proceeding with the massive gravity theory, in- 
dices are raised/lowered, from now on, with the dynam- 
ical metric g^i/; for example = g^^''Hpv. 
We define a new tensor /C",, as 

^,'' = '^/-(v/.9-M.9-i?]) \ (4) 

where the square root is formally understood as 
vG4^ y/A^ = . This allows us to represent the com- 
plete potential for gravitational interactions as 

C = ^^g{R- m'U) , U=[tv (^2) _ (tr ]Cf] . 

(5) 

By expanding the potential in H^^, we get an infinite 
sum of interaction terms for H^^, with the FP term at 
the lowest order. 

In order to study exact solutions associated with the 
previous Lagrangian, it is convenient to express /C in 
terms of matrices, namely 

= I - ^g-^v+Ql (6) 

where I denotes the identity matrix, and we have used 
H^v = QtLv - iVfj-iy + Qfj.iy)- The potential in four dimen- 
sions then reads 

U tT g-^ [r] + Q] - 12 

+ tr [ri + Q] (6 - tr Vg-'iv+Q]) ■ (7) 

In general, the task is to calculate the trace of 
\/g^^r]+Q\. Given that g~^ [rj + Q] is a square ma- 
trix, the Schur decomposition theorem ensures that it 
can be expressed as 

g-'[v + Q] =rDr-\ (8) 

for some unitary matrix T, and an upper triangular ma- 
trix D. The diagonal entries of D are the eigenvalues 
of g~^ [r] + Q], and we call these eigenvalues Ai, . . . , A4. 
Then, since g~^ [i] + Q] — T ^pDT^^ ^ one can express 
the traces in the formulae above, in terms of eigenvalues, 
as 

trg-M^Z+Q] = 

i 

trv/5-M^/+S] = Ev^- 



Plugging these expressions into the potential, Eq. (7), 
we find the following expression for U 

" = E^^ + (Ev^)(6-Ev^) -12- (10) 

i j i 

Spherically symmetric configurations: We now 

focus on analyzing properties of spherically symmetric 
configurations in this set-up. We start by considering 
static configurations in the unitary gauge, tt^ = (see 
Ref. [15] for spherical symmetric solutions in the FP the- 
ory). The most general form of the metric respecting 
spherical symmetry is 

ds^ = -C(r) dt^ + A{r) dr'^ + 2D{r) dtdr + B{r)dn^, 

(11) 

where dO^ = dO-^ + sin^ 9d(fP . We choose to write the 
non-dynamical flat metric as ds^ = —dt^ + dr^ + r^d^l'^. 
Notice that in GR one can set B{r) = r-^ by a coordi- 
nate transformation, but this is not possible here, since 
we have already fixed the gauge completely. In order to 
simplify the analysis, it is convenient to define the combi- 
nation A (r) = A{r)C{r) + D'^{r). We plug the previous 
metric into the Einstein equations 

Gf^u = T^„, (12) 

where the energy momentum tensor from the potential U 

ofEq. (10) is defined as T^^ = The Einstein 

tensor G^u satisfies the identity D{r) Gtt + C{r) Gtr = 0, 
which implies the algebraic constraint 

= D{r)TY,+C{r)TY,. 

D{r) (2r - 3^B{^) 

= 7^ 7Tr2- (13) 

[A{r)+C{r)+2^/W)) 

The previous condition can be satisfied in two ways, 
which lead to two different branches of solutions. We 
can either set D{r) ~ 0, and focus on diagonal metrics, 
or alternatively, set B{r) = 4r^/9. The fact that there 
are two branches of solutions indicate that, unlike in GR 
where Birkhoff's theorem holds, there is no uniqueness 
theorem for spherically symmetric solutions in this the- 
ory. We analyze the two branches in turn. 

Diagonal-metric branch - asymptotically flat so- 
lutions: The case of diagonal metrics, D{r) = in Eq. 
(11), leads to equations which in general cannot be solved 
analytically. Thus, we will analyze them perturbatively, 
showing that they lead to asymptotically flat solutions 
(see Ref. [16] for a more detailed discussion on this branch 
of solutions). We expand the functions A, C and B as 

A{r) - C{r) = {1+nf, B{r) = j^^^, (14) 

and truncate the field equations to first order in n, / 
and h. It is more convenient to introduce a new radial 
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coordinate p = B{r), so that the hnearised metric is 
expressed as 

ds^ = -(1 + 2n)dt^ + (1 - f)dp^ + p^dn^, (15) 

where / = / — 2/i — 2ph' and a prime denotes a derivative 
with respect to p. In this new coordinate, the solutions 
for n and / are then given by 



2n 



8GM 
AGM 



(1 + mp)e 



(16) 



where we fix the integration constant so that M is the 
mass of a point particle at the origin and SttG = M^^^. 
See the right plot of Fig. 1 for the general behaviour of 
these solutions. Notice that these configurations, as an- 
ticipated, are asymptotically fiat and exhibit the vDVZ 
discontinuity, i.e. they do not agree with the GR solu- 
tions {2ncR = Jgr = —2GM/p) in the limit m ^ 0. 
However, in order to understand what really happens in 
this limit, one should take into account the non-linear be- 
haviour of h. Let us consider scales below the Compton 
wavelength mp ^ 1, and at the same time ignore higher 
order terms in GM. Under these approximations, the 
equations of motion can still be truncated to linear order 
in / and n, but since h is not necessarily small, we keep 
all non-linear terms in h. Then we obtain the following 
equations 



2pn' = 

/ = 
GM 



2GM 



{mpfh, 



P 

-2^-{m^pf{h-h% 



(17) 



We should stress that these are exact equations in the 
limit mp <^ 1, GM/ p <^ 1, i.e. there are no higher or- 
der corrections in h. For large radial p values, one can 
linearise the equations in h, recovering the solution in 
Eqs. (16), to first order in mp. On the other hand, the 
Vainshtein mechanism applies, and below the so-called 
Vainshtein radius, pv = (GMm~^)^/'^. h becomes larger 
than one and the non-linear terms in h become impor- 
tant, recovering GR close to a matter source. Actu- 
ally, for p ^ Pv the solution for h is simply given by 
\h\ = pv/p^ 1- The latter solution for h and Eq. (17) 
imply 



2pn' 
f 




(18) 



Therefore, corrections to the GR solutions are indeed 
small for p < py, as shown in the left plot of Fig. 1. 



The Vainshtein mechanism becomes also transparent 
in a non-unitary gauge. Indeed by performing the coor- 
dinate transformation p — ^ B{r), we excite the p com- 
ponent of the Stiickelberg field (see Eq. (3)), ttP = —ph. 
Thus the strong coupling nature of h is encoded in ir^ 
in this coordinate. It is possible to construct an effec- 
tive theory for this Stiickelberg field in the so-called de- 
coupling limit [13]: first we introduce a scalar so that 
T^fj. = dfj,TT/Al, where A3 = m?Mpi. Then the covarianti- 
zation of metric perturbations H^y is written as 



Mpim^ 



(19) 



IT FT" 



Formally, 



where 11^^ = 9^a^7r and 11^^, 
the decoupling limit is achieved by taking m — > and 
Mpi —J- 00, but keeping A3 fixed. By substituting Eq. (19) 
into the action, one can show that the kinetic terms of tt 
become total derivatives and a mixing appears between 
/i^^ and TT, which can be diagonalised using the definition 



(20) 



The Lagrangian is then written as 

C = £GR(V) + ^^n7r-^(9^)2n^ 

+ ^(57r)^([n']-[n]') 

+ ^(a7r)2([n]3-3[n][n2] + 2[n3]), (21) 

where [H] = H^j, [H^] = n^-^n^,, [H^] = n'-'^n.^n^ 

and £gr is the linearised Einstein- Hilbert action for h^^. 
The terms containing the scalar field tt are known as 
Galileon terms, which give rise to the second order dif- 
ferential equations, as explained in the Introduction. For 
the spherically symmetric case, the equation of motion 
for TT simplifies to [12, 17] 
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Af 



M 



4tt Mpi p3 



(22) 



where the integration constant is again chosen so that M 
is a mass of a particle at the origin. Using the relation 
between tt and h, h = — tt'/ (m^Mp/p), we can show that 
the solutions for /, n and h given by Eq. (17) agree with 
the solutions Eq. (20) and Eq. (22). 

We have shown that the weak field solutions for the 
metric Eq. (11) with D{r) = have three phases (see 
Fig. 1). On the largest scales, ^ p, beyond 

the Compton wavelength, the gravitational interactions 
decay exponentially due to the mass of graviton: see 
Eq. (16) and region 3 in Fig. 1. In the intermediate region 
Pv < p < m^^, we obtain the 1/r gravitational potential 
but Newton's constant is rescaled G — > 4G/3. More- 
over, the post-Newtonian parameter 7 is 7 = f /{2n) = 



4 



(1/2)(1 + mp), which reduces to 7 = 1/2 in the m — )■ 
limit, instead of 7 1 of GR, showing the vDVZ dis- 
continuity (see region 2 in Fig. 1). Finally, below the 
Vainshtein radius p < pv, the GR solution is recovered 
due to the strong coupling of the tt mode (see Eq. (18) 
and region 1 in Fig. 1). This background solution pro- 
vides us with a testing ground for the BD ghost. Instead 
of expanding the action in H/^^i, around the Minkowski 
spacetime perturbatively, one can study linear pertur- 
bations around this non-perturbative solution using the 
complete potential Eq. (7). In order to obtain the fully 
non-linear solution, a numerical approach is necessary. 
In the next section, we consider the second branch of 
solutions for this theory, which can instead be obtained 
analytically. 




FIG. 1: Numerical solution for drf = /', drfi — n' , and the 
quotient 7' = f /2n around the Vainshtein radius pv (left) 
and the Compton Wavelength p ~ 1/m (right). Region 1 (2) 
shows how GR solutions are (not) recovered inside (outside) 
the Vainshtein radius pv- Region 3 shows the asymptotic 
decay of the linear solutions (Eq. (16)). Here, GM — 1. 

Non-diagonal-metric branch - Schwarzschild de 
Sitter solutions: Next, we analyze the second branch 
of vacuum solutions that solve Eq. (13), where B{r) = 
4r^/9. Interestingly, this branch leads to asymptoti- 
cally de Sitter configurations. There is another identity 
C{r)Tl^^ + A{r)TYt = 0, which leads to the condition 

A(r) = A{r)C{r) + D^{r) = Aq = const. (23) 

The remaining Einstein equations provide the following 
unique solution (see Ref. [16] for detailed derivations) 

A{r) = ^{p{r)+a + l), B{r) = (24) 



C{r) 
where 



9Ao 



9 



^ il-p{r)), D{r) = -AoVpWbW+«), 



2p rn^r'^ 16 



1, 



(25) 



with arbitrary p and Aq. This solution is similar to that 
in [18], up to numerical factors. Notice that this configu- 
ration depends on two integration constants. A sufficient 



condition to ensure that D{r) is real, is to choose /i > 
and < \/Ao < 4/9. The form of metric coefficients as in 
Eq.(24) do not allow a manifest comparison with de Sit- 
ter spacetime, since we have already chosen the unitary 
gauge and cannot do a further coordinate transformation 
without exciting components of tt^. However, if we al- 
low for a vector tt^ of the form tt^ = (7ro(r), 0, 0, 0), the 
metric can be rewritten in a diagonal form as 



-C{r) + A{r) dr^ -t- B{r)dVt^ 



(26) 



Then we can write down the action in terms of C, A, B 
and ttq, considering them as fields. It is possible to show 
that the following configuration solves the corresponding 
equations of motion 



= n 



9 l-p{r)' 



while C(r) and B{r) are the same as in Eq. (24). 
The resulting metric has then a manifestly de Sitter- 
Schwarzschild form by making a time rescaling t — > 

1 /2 

(4/9Ao' )t. However, we should note that this time 
rescaling cannot be done without introducing an ad- 
ditional time dependent contribution to ttq. As ex- 
pected, the metric in Eq. (26) can be obtained by mak- 
ing the following transformation of the time coordinate 
dt = dt + n'^dr to the metric (11); this produces a non- 
zero time component of tt'^, that does not vanish even in 
the TO limit for any allowed value of Aq. There are 
two integration constants, p and Aq, in this solution. In 
GR; p corresponds to the mass of a source at the origin 
but a careful analysis including a matter source is nec- 
essary to fully understand the role of these integration 
constants. Note that there is an apparent singularity at 
the horizon p(r) = 1, both for the metric and ttq. 

We can further make a coordinate transformation at 
the expense of exciting further components of tt'^. For 
example, by setting pi = and making the following co- 
ordinate transformations t = Ft{T, p), r = F^[t, p) with 



Ftir,p) 



Fr{T,p) 



4 /sinh(2^) -f ^;V^e"^/2 
arctanh ' ^ ^ ' 



3Ay^m 
3 



Vcosh(i^) - IIiiele™-/2 ) 

mr/2 ^28) 



the metric becomes that of flat slicing of de Sitter, 

ds2 ^ -dT^ + e'^^dp^ + p^dn'^), (29) 

where the Hubble parameter is given by m/2. The 
Stiickelberg fields tt^ are now given by tt^ = 

iT:^T,p),TTP{T,p),0,0), TT- = + T - Ft (t, p) , Tt" = 

p — Fr (r, p) . This is an interesting solution in which 
the acceleration of the universe is determined by the 
graviton mass and the Hubble parameter is given by 
to/2. For Aq = 16/81, this solution reduces to the "self- 
accelerating" solution obtained in the decoupling limit in 
Ref. [19]. 
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Conclusions: The solutions obtained in the non- 
linear covariant massive gravity are remarkably similar 
to those in the DGP braneworld model including the ex- 
istence of the "self-accelerating" de Sitter solution with- 
out cosmological constant [20] although there are differ- 
ences in detail. There are a number of important issues. 
Firstly, we should confirm that there is no BD ghost in 
this theory by studying perturbations around the non- 
perturbative solution obtained in this letter. In the DGP 
model, the self-accelerating solution suffers from a ghost 
instability [10, 11, 21], which is related to the ghost in 
the FP theory on a de Sitter background. It is crucial to 
study the stability of the de Sitter solution in this model. 
In fact Ref. [19] showed that there exists a ghost in this 
self-accelerating background in the decoupling limit for 
a particular value of the second integration constant Aq. 
They argue that this ghost can be cured by adding higher 
order corrections in /C to the potential. Our formalism is 



ready to be applied to this extended model. However, we 
believe a more complete analysis of perturbations about 
our exact solution is needed. Once these issues are clar- 
ified, the massive gravity model presented here provides 
an interesting playground to study large distance modi- 
fications of general relativity. 
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